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Abstract 

We study a monomer-dimer model with repulsive interactions between the 
same species in one dimension. With infinitely strong interactions the model 
exhibits a continuous transition from a reactive phase to an inactive phase 
with two equivalent absorbing states. Monte Carlo simulations show that 
the critical behavior is different from the conventional directed percolation 
universality class but seems to be consistent with that of the models with the 
mass conservation of modulo 2. 

PACS numbers: 64.60.-i, 82.65. Jv, 02.50.-r, 05.70.Ln 
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A monomer-dimer model was introduced by Ziff, Gulari, and Barshad to describe the 
oxidation of carbon monoxide on catalytic surface [Q. In this model, a monomer (CO) 
adsorbs onto a single vacant site, while a dimer (O2) adsorbs onto a pair of adjacent vacant 
sites and then immediately dissociates. A nearest neighbor of adsorbates, comprised of a 
dissociated O atom and an CO atom, reacts and forms a CO2 molecule and desorbs from the 
metal surface. In two dimensions, as the CO gas pressure is lowered, the system undergoes 
a first-order transition from a CO-saturated inactive phase into a reactive steady state and 
then a continuous transition into a 02-saturated inactive phase. This continuous transition is 
shown to be in the same universality class as the directed percolation (DP) and the Reggeon 
field theory (RFT). 

Motivated by the monomer-dimer model, there have appeared many related lattice mod- 
els to study nonequilibrium phase transitions. One of the simplest models is the monomer- 
monomer model in which particles of two different species can adsorb on a single vacant 
site. If two adsorbates of different species are nearest neighbors each other, they react and 
form a product which desorbs immediately from the surface. The system exhibits a first- 
order transition from a phase saturated with one species to another. Allowing desorptions 
of one species leads to a continuous transition which belongs to the DP universality class 




In both the monomer-monomer and the monomer-dimer models, interactions between 
adsorbing particles are ignored except for the actual surface reaction which converts the 
reactants into the product. Recently, an interacting monomer-monomer model p,^] has been 
introduced where particles of the same species have variable repulsive interactions. When 
the interaction strength is weak, the interaction model exhibits only a first-order phase 
transition between two saturated phases. At the critical value of the interaction strength, a 
first-order line terminates at a tricritical point beyond which two continuous-transition lines 
appear. The two saturated phases are separated from a reactive steady state by these two 
lines. These two continuous transitions are shown to be in the DP universality class again. 

A common feature of these models is that they exhibit, if any, a continuous transition 
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from a reactive phase into an inactive phase of a single absorbing state. The resulting 
critical behaviors are classified into the category of the so-called "DP conjecture" 



which depicts that models exhibiting a continuous-transition to a single absorbing state 
generally belong to the universality class of the directed percolation. 

The universality class for models with a single absorbing state is well established. But 
few studies have been made for models with more than one absorbing states. Very recently 
Jensen and Dickman ||ll],|l2| have extensively studied some nonequilibrium lattice models 



with infinitely many absorbing states, the pair contact process (PCP) and the reaction dimer 
(RD) model. Both models have a continuous transition from a reactive phase into an inactive 
phase with infinitely many absorbing states, which is shown to be, rather surprisingly, again 
in the DP universality class. 

In contrast to the well-established DP universality class, only a few models have appeared 
in literature which are known to be in the different universality class from DP. Grassberger, 



Krause, and von der Twer p^JT5[] studied two models of probabilistic cellular automata 
(PGA), namely model A and B. Both models exhibit a continuous transition from an active 
phase into an absorbing phase, in which the system is trapped into one of two (translation- 
ally) symmetric states with particles and vacancies placed in the alternating sites. But these 
models behave differently in the absorbing phase. Once the system enters into one of the 
two absorbing states, it remains in that state forever in model A, but oscillates from one 
state to the other in model B. In spite of the discriminating behavior in the absorbing phase, 
both models are shown to be in the same universality class but different from DP. The order- 
parameter exponent (3 was obtained by static Monte Garlo simulations, (3 = 0.6 ± 0.2 ||14|| , 
but later by dynamic Monte Garlo simulations, (3 = 0.93 ± 0.12 utilizing the hyperscaling 
relation The evolution rules for both models in common involve the processes X —>■ 3X 
and 2X — > that the number of particles is conserved modulo 2. Grassberger has addressed 
in [|I5| that this conservation law might be responsible for the non-DP behavior. 

Another model which conserves the number of particles by modulo 2 can be found in 
branching annihilating random walks with even number of offsprings (BAW). Even though 



the BAW has a single absorbing state, static Monte Carlo simulations [|T6| show clearly 
non-DP behavior with {3 = 0.7 ± 0.1 in one dimension. This value seems to be consistent 
with that of the PCA, even though heavy numerical errors obscure the fact. The critical 
behavior of BAW with odd number of offsprings is compatible with DP . 

All previous works |T^-|TB|,|TT],|T^ seem to imply as a whole that more important in de- 
termining the universality class of nonequilibrium phase transitions is not the number of 
absorbing states but the mass conservation law in dynamics. This shows a sharp contrast 
to the case of equilibrium phase transitions where the number of ground states is relevant 
to the universality class. The symmetry between absorbing states may be more important 
than the number of absorbing states in determining the universality class. In the PCP, the 
infinitely many absorbing states are not equivalent probabilistically. Some absorbing states 
can be reached more easily than other absorbing states by the PCP dynamics. Therefore it 
is important to study a model with multiple equivalent absorbing states. The PCA models 
have two equivalent absorbing states but also another conserving quantity, the mass con- 
servation of modulo 2, which is claimed to be responsible for the non-DP behavior. In this 
letter, we address a question whether a model with multiple equivalent absorbing states 
may have non-DP behavior without any mass conservation law. As one of the simplest such 
models, we study the monomer-dimer model with infinitely strong repulsive interactions 
between the same species. 

The interacting monomer-dimer model (IMD) is a generalization of the simple monomer- 
dimer model, in which particles of the same species have nearest-neighbor repulsive inter- 
actions. This is parameterized by specifying that a monomer (A) can adsorb at a nearest- 
neighbor site of an already-adsorbed monomer (restricted vacancy) at a rate r^fc^ with 
< ta < 1, where is an adsorption rate of a monomer at a free vacant site with no 
adjacent monomer-occupied sites. Similarly, a dimer {B2) can adsorb at a pair of restricted 
vacancies {B in nearest-neighbor sites) at a rate r^/cs with < < 1, where is an ad- 
sorption rate of a dimer at a pair of free vacancies. There are no nearest-neighbor restrictions 
in adsorbing particles of different species. The case = = 1 corresponds to the ordinary 



noninteracting monomer- dimer model which exhibits a first-order phase transition between 
two saturated phases in one dimension. In the other hmiting case = = 0, there exists 
no fully saturated phase of monomers or dimers. However, this does not mean that this 
model has no absorbing states any more. In fact, there are two equivalent absorbing states 
in this model. These states comprise of only the monomers at the odd- or even-numbered 
lattice sites. There needs a pair of adjacent vacancies for a dimer to adsorb, so a state with 
alternating sites occupied by monomers can be identified with an absorbing state. 

In this letter, we consider the IMD with = = only for simplicity. General 



phase diagram of the IMD will appear elsewhere . Then the system can be characterized 
by one parameter p = ky^/{kA + ks) of the monomer adsorption-attempt probability. The 
dimer adsorption-attempt probability is given by g = 1 — p. The order parameter of the 
system is the concentration of dimers p in the steady state, which vanishes algebraically as 
p approaches the critical probability pc from below: 

P ~ (Pc - pf , (1) 

where f3 is the order-parameter exponent. There are a characteristic length scale ^ and time 
scale r which diverge at criticality as 

e~be-pr^S (2) 

r~be-p|~^ii, (3) 

where z/^ (z/||) is a correlation length exponent in the space (time) direction. It is quite 
difficult to measure the order parameter accurately near criticality by static Monte Carlo 
simulations due to the critical slowing down. Moreover, there are strong finite-size effects 
near criticality because of the diverging correlation length. In this letter, we utilize the 
finite-size scaling (FSS) idea developed for nonequilibrium phase transitions by Aukrust, 
Browne, and Webman pO[1 . 

Various ensemble-averaged quantities depend on system size through the ratio of the 
system size and correlation length L/^. Thus we can take the concentration of dimers near 
criticality as the following form: 



p(p,L)r.L-^/^-f((p,-p)Ly^-) (4) 

such that at Pc 

p(p,, L) ~ L-f''^- (5) 

and 

f{x) as X — > oo. (6) 

In the supercritical region (p < pc), the concentration p remains finite in the hmit L — > oo, 
but it should vanish exponentially in the subcritical region (p > pc) ■ 

For the time dependence of the order parameter at criticality, one may assume a scaling 
form 

-p{L,t)^L-PI''-g{Lt-''-'''\\), (7) 

so that for short time (or L — > oo) p has a power-law dependence on t as 

pr^r'^'^'K (8) 

The characteristic time r for a finite system is defined as the elapsing time for a finite system 
to enter into the (quasi-) steady state. An elementary FSS analysis finds r at Pc 

T{p,,L)^m'^\ (9) 

We run static Monte Carlo simulations. The initial configuration is far from absorbing 
states with all sites vacated and we use periodic boundary conditions. Then the system 
evolves along the dynamical rules of our model. After one adsorption attempt on the average 
per lattice site (one Monte Carlo step), the time is incremented by one unit. The system 
reaches a quasi-steady state first and stays for a reasonably long time before finally entering 
into an absorbing state. We measure the concentration of dimers in the quasi-steady state 
and average over some independent samples which have not yet entered the absorbing states. 
The number of independent samples varies from 50,000 for system size L = 32 to 300 for 
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L = 2, 048. The number of time steps ranges from 500 to 2 x 10^ and, at least, 200 samples 
survive until the end of simulations. 

At pc, we expect the ratio of the concentrations of dimers for two successive system 
sizes p(L/2)/p(L) = 2^^'^^, ignoring corrections to scaling. This ratio converges to unity for 
p < Pc and diverges to infinity for p > pc in the limit L — oo. We plot the logarithm of 
this ratio divided by log 2 as a function of p hi L = 64, 128, 256, 512, 1,024, and 2,048 
in Fig.l. This plot shows strong corrections to scaling. The crossing points between lines 
for two successive sizes move to the right as the system size grows. In the limit L ^ oo, 
we estimate the crossing points converge to the point at pc — 0.581(2) and /3/z/_l ^ 0.45(2). 
The critical probability Pc can be more accurately estimated from dynamic Monte Carlo 
simulations |jT9|. The value of the exponent ratio f3/i'± is clearly different from the standard 
DP value of 0.2524(5) which is combined with z/^ = 1.0972(6) and (3 = 0.2769(2). 

In Fig. 2, we have plotted pL^/'^^ against x = {l—p/pc)L^^'^^ in a double-logarithmic plot. 
From Eqs.(^ and (^, it follows that for small x the curve converges to a constant value, 
while for large x it should follow a straight line with slope (3. We find that the data for the 
various system sizes are best collapsed to a curve with choices Pc = 0.581(1), = 1.73(3), 
and (3/i'± = 0.45(1), from which we get /3 = 0.78(3). The straight line in Fig.2 is a line with 
slope 0.78, which is in excellent accord with asymptotic behavior of x in Eq.(^. 

By analyzing the decay characteristic of dimer concentrations p{t,L), we can determine 
the value of another correlation exponent z/y . First the finite-size behavior of the character- 
istic time T is investigated. At p = Pc — 0.581, we measure r for the various system sizes. 
The double-logarithmic plot (see Fig. 3) for the characteristic time r versus the system size L 
shows a straight line from which we obtain the value z^||/z^± = 1.664(3). This gives the value 
of z/|| = 2.88(5), combined with the value z/_|_ = 1.73(3). We also can check the value of z/y 
independently by measuring the decay exponent of p at p = Pc- From the double-logarithmic 
plot for p(t) versus the time t in Fig. 4, we estimate P/i'w = 0.276(2), which is in excellent 
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agreement with the above results. We summarize our results for the critical exponents as 

(5 = 0.78(3), = 1.73(3), i/y = 2.88(5). 

In summary, we have numerically studied the interacting monomer-dimer model with 
infinitely strong repulsive nearest-neighbor interactions between the same species in one 
dimension. The system exhibits a continuous transition from a reactive phase into an inactive 
phase with two equivalent absorbing states. This model does not involve the dynamical 
processes with the mass conservation of modulo 2. Nevertheless, the values of the critical 
exponents are clearly different from the DP values and seems to be consistent with the 
values of the PCA and the BAW, in which the the mass conservation of modulo 2 governs 
the dynamics. Therefore we conclude that the symmetry between absorbing states as well 
as the mass conservation law are equally important in determining the universality class for 
nonequilibrium phase transitions. Up to date, it is not yet clear why the model with two 
equivalent absorbing states and the one with the mass conservation of modulo 2 seem to be 
in the same universality class. Dynamical exponents and hyperscaling relations of the IMD 
will be discussed elsewhere |TB| . 
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Figure Captions 

Fig. 1 : Plots of log [p(L/2)/p(L)] / log2 versus p. Open squares are for L = 64, solid 
squares for L = 128, open circles for L = 256, closed circles for L = 512, open 
triangles for L = 1, 024, and closed triangles for L = 2, 048. 

Fig. 2 : The double-logarithmic plot for the data of pL^/'^^ against a; = (1 —p/pc)L^^'^^ for 
the various system size L = 32 - 1, 024. With p^ = 0.581, z/_l = 1.73, = 0.45, the 
data are collapsed to a curve. The solid line is of slope 0.79 (= P). 

Fig. 3 : The characteristic time against the system size L in a double-logarithmic plots for 
the various system size L = 32 — 1, 024. The solid line is of slope 1.664 (= z/y/z/^). 

Fig. 4 : The time-dependence of the concentration of dimers at pc = 0.581. The slope of 
the curve gives the value of = 0.276(2). 
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